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Abstract. We have measured p t -dependent two-particle number correlations on azimuth 
and pseudorapidity for eleven centralities of ^/snn = 62 and 200 GeV Au-Au collisions at 
STAR. 2D fits to these angular correlations isolate the azimuth quadrupole amplitude, denoted 
2v2{2D}(pt), from localized same-side correlations. Event-plane V2(pt) measurements within the 
^ ' STAR TPC acceptance can be expressed as a sum of the azimuth quadrupole and the quadrupole 

D . component of the same-side peak. V2{2D}(p t ) can be transformed to reveal quadrupole pt 

spectra which are approximately described by a fixed transverse boost and universal Levy form 
nearly independent of centrality. A parametrization of V2{2D}(pt) can be factored into centrality 
and pt-dependent pieces with a simple pt dependence above 0.75 GeV/c. Results from STAR 
are compared to published data and model predictions. 



=5 



> 

m 

^ ■ 1. Introduction 
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^-j. , The collective behavior of particles produced in heavy-ion collisions is one of the major topics 

of study at RHIC, and a large quadrupole component in the distribution on azimuthal angle 
<j) has been observed [TJ. The large quadrupole azimuthal distribution in non-central collisions 
is usually described by elliptic flow, in which the different pressure gradients in and out of 
plane turn the initial position-space anisotropy of the colliding nuclei into a momentum-space 
anisotropy [2J. 

The azimuthal anisotropy is typically measured with the quantity V2, the second Fourier 
component of the distribution of particle emission on azimuth with respect to the angle of the 



C3 



reaction plane [3]. In real events this reaction plane angle is not known so many methods exist 
to estimate V2- The present analysis uses a method based on fitting two-dimensional angular 
correlations similar to [I] but expanded to include the dependence on transverse momentum as 
well as centrality. A simultaneous analysis of centrality and transverse momentum dependence 
of V2 in the absence of "nonflow" from the same-side 2D peak (jet-like correlations) is presented. 

One method of measuring v% is the event-plane method, ^{EP}, in which particles in the 
region of interest are used to estimate the reaction plane [3]. A related approach, t>2{2}, 
calculates V2 directly from two-particle correlations [5]. These two methods are essentially 
equivalent (^{EP} ~ ^{2} to 5%) [6]. This analysis, denoted by «2{2D}, uses an approach 
similar to t>2{2} but includes the dependence of the two-particle distribution on pseudorapidity 
as well as azimuth. Significant jet-like correlations are observed at all centralities [7J [8] and 
their 2D angular dependence can be used to distinguish them from the azimuth quadrupole 
component, which within the STAR TPC acceptance has no significant r/ dependence [9]. Jet- 
like correlations are separated from ^2{2D} by 2D model fits described in Sec. [3j 



Figure 1. Left Panel: Two examples of marginal distribution cuts in yt x yt space. Remaining 
panels: An example of the yt evolution of correlation structures for 40-50% central collisions. 
The plots correspond to yt bins of 1.4 < y t < 1.8, 3.0 < yt < 3.4, and 3.8 < yt < 4.2. 

2. Correlations 

Angular correlations without a trigger particle are constructed by considering all pairs of particles 
in an event. For each particle the azimuth angle (ft, pseudorapidity r/ and transverse momentum 
Pt are measured, defining a six dimensional two-particle space. Integrating over pt and projecting 
the angular variables onto their difference axes t/a = Vi ~ V2 and (ft a = 4>i ~ 4*2 yields 2D angular 
autocorrelations [H 13 [TU] . 

The correlation measure is constructed from pair densities. Sibling pairs from the same 
event contain the primary signal. Mixed pairs from different events are useful as a reference 
for removing detector artifacts. The basic correlation measure we use is Ap/p re f, the ratio of 
the sibling minus the mixed pairs to the mixed pairs, a per-pair measure. It is also sometimes 
desirable to use a two-particle measure like Ap/^/p~^J [11 1. However, this analysis uses Ap/p re f 
which is easily converted to v 2 '- 1v\\2D} = Ap[2]/ p re f, where the [2] indicates the quadrupole 
term of the fit model in Eq. [T] below. 

When we consider momentum dependence it is desirable to use transverse rapidity instead 
of transverse momentum: yt = In {(pt + mtj/mo}. Transverse rapidity is the analogue to 
longitudinal rapidity but in the transverse direction. For unidentified particles we make the 
approximation mo = to define a regularized ln(pt) measure. 

We construct marginal distributions by restricting the yt range of one of the particles in the 
two-particle correlation to a certain interval. Due to the inherent diagonal symmetry in yt x yt 
space this corresponds to making cross-shaped cuts as seen in Fig. [TJ (left panel). These cuts are 
analogous to the procedure in a standard event-plane ^(Pt) analysis [3] in which particles from 
the full momentum range are used to estimate the event plane and particles in the momentum 
interval are correlated to it. The marginal distributions allow us to study correlations at large 
Pt values and are easily converted into the usual v^ipt) measure for comparison with published 
v 2 {pt) data. 

In this analysis there is a minimum pt cut for all particles of 150 MeV/c which corresponds to 
a pion yt of 0.93. There are nine yt bins. The first yt bin includes all particles with yt less than 
1.4 and the seven subsequent bins are 0.4 units of yt wide. The final bin includes all particles 
with y t greater than 4.2. This corresponds to particles above ~ 5 GeV/c. 

3. Fit Parameters 

This analysis is based on 14.5 million Au-Au collisions at y / sjv r /v = 200 GeV and 6.7 million 
Au-Au collisions at ^/snn = 62 GeV observed with the STAR TPC. The acceptance was 
defined by transverse momentum pt > 0.15 GeV/c, \r]\ < 1 and 2ir azimuth. Minimum-bias 
event samples were divided into 11 centrality bins: nine ~ 10% bins from near 100% to 10%, 
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Figure 2. The amplitude, width in ?7a, and width in 0a of the 2D Gaussian and the amplitude 
of the quadrupole component of the fit for Au-Au 200 GeV 60-70%, 50-60%, 40-50%, 30-40%, 
20-30%, and 10-20% collisions. Error bars are for fit errors only. 



the last 10% divided into two ~ 5% bins. 

Two-particle correlation histograms are constructed for each of the 11 centrality and 9 yt bins 
described above for a total of 99 bins at each ^/snn- Example histograms are shown in Fig. [TJ 
Fig. [2] shows the evolution of the most relevant fit parameters at ^/saT/v = 200 GeV with yt for 
six centrality bins. Systematics errors are still under study. 62 GeV results have similar trends 
as 200 GeV. 

In order to distinguish the quadrupole signal from ^-localized (jet-like) correlations we will 
fit the major structures in the 2D histograms. The fitting model used is similar to that used 
in [1], but the sharp peak from electron pairs is excluded rather than fit. The histograms are 
constructed with 24 bins on 0a with centers ranging from —ir/2 to 177T/T2 in steps of 7r/12 and 
25 bins on t?a with centers ranging from —1.92 to 1.92. A total of seven bins are excluded from 
the fit, where f?A = and 0a = 0, ±7r/12 and where 0a = and t/a = ±0.08, ±0.16. 

The complete fitting function is then 

F = ^ A cos(0A)±^ A cos(20 A )±^oe + A x e l V ^ y V vaJ J + A 2 . (1) 

There are 5 terms and 8 parameters. Two of the terms are sinusoids on 0a , one is a constant 
offset, and one models only structure on r/A- Thus, the remaining term — the two-dimensional 
Gaussian — is solely responsible for describing "nonflow" [6]. 

Fig. E] (first three panels) shows measured quadrupole results from this analysis vs. pt 
compared to previously published ^{EP} results for 30-40%, 5-10% and 0-5% central collisions. 
In the case of 0-5% the quadrupole component was found to be too small to measure reliably so 
the cross-hatched region represents an upper bound. We calculate the contribution of the 2D 
Gaussian peak to the second Fourier component using the measured amplitudes and widths in 
the curve labeled 'jets'. We find V2{EP} seems to be dominated by the SS peak structure in 
5-10% collisions and entirely due to it in 0-5% collisions. This has important implications for 
the interpretation of published V2 data. 

4. Quadrupole Descriptions 

Consider a model in which the single-particle density on yt and and corresponding transverse 
boosts can be decomposed into azimuth-independent po and quadrupole P2 pieces [12]. Define 
V2 (j/t , 6) = po(yt,b)v2(yt,b) and invoke a blastwave model to get the relation V2{yt,b) « 
PtAyt2(b)p2(yt, b)/{2T2) where Ay t 2 is the quadrupole component of the radial boost distribution, 
P2 is the quadrupole spectrum, and T2 is a blastwave parameter. Now construct the unit-integral 
quantity Q(y t ,b) = [V 2 {yt, b)/pt]/[V2(b)(l/p t )} « [p2(yt, b)}/[p 2 (b)}. Most of the parameters drop 
out and this variable directly relates measurable quantities to the quadrupole spectrum model. 
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Figure 3. First three panels: Comparison of published ^{EP} for Au-Au 200 GeV 30-40%, 
5-10% and 0-5% central collisions [9j and results from this analysis. Errors are fit errors only. 
Right Panel: The quadrupole spectrum data Q(yt,b) and fit Qo(yt)- The dotted curve is a 
zero- viscosity limit of a viscous hydro theory defined in |13j . 



There is little centrality dependence to the quantity Q(yt,b). A single Levy distribution 
Qo(yt) can describe the data within 10% with a few simple parameters: a yt boost Ayto = 0.58 
and parameters T2 = 0.09 GeV and n<i = 13.8. There is deviation from this curve for very 
central collisions at low-p^. A similar analysis was made in [12] for minimum-bias identified 
particles without centrality dependence and the observed boost values agree. It is notable that 
the boost is a single sharp value and not a distribution as would be expected from transverse 
Hubble expansion. 

A related and more accurate parametrization has been made. At higher pt it is observed that 
V2{2D}(pt , b) is proportional to ptexp(— pt/4) so we derive the form: 

v 2 {2B}(p t ,b) « (l/pt)^{2D}(6)p t exp(-p t /4) x f( Pt ,b), (2) 

where {1/pt) ~ 2.1 (GeV/c) -1 and f(pt, b) is a dimensionless factor necessary for the full pt and 
centrality dependence. It is fit to data to get f(pt,b) = 1 + C(6)[erf(y 4 — 1.2) — erf (1.8 — 1.2)] 
and C(b] ) = 0.12 - {v - 3.4)/5 - [{v - 3.4)/2] 5 . 

Using 2D correlations we observe V2(pt) for a wide range of centralities without significant 
"nonflow" contamination and conclude that published V2(pt,b){EP} results [9] contain a 
significant contribution from the same-side peak correlations. 

We have measured the quadrupole on y t as a function of b and deduced the quadrupole 
spectrum variation with centrality. This is the first complete factorization of the scaling 
properties of the quadrupole. A detailed study of the systematics is in progress. When the 
quadrupole spectrum is isolated we also observe a single 6-independent yt boost value that is 
consistent with the minimum-bias results for different particle species |12j . We developed an 
accurate parametriazation of -U2{2D} over a large kinetic domain and surprisingly simple trends 
are observed. 
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